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Abstract. A formal fourth order differential operator with a singular coefficient 
that is a linear combination of the Dirac delta-function and its derivatives is con- 
sidered. The asymptotic behavior of spectra and eigenfunctions of a family of 
differential operators with smooth coefficients approximating the singular coeffi- 
cients is studied. We explore how behavior of eigenvalues and eigenfunctions is 
influenced by singular coefficients. The limit operator is constructed and is shown 
to depend on a type of approximation of singular coefficients. 



1. Introduction 

Differential operators with singular coefficients appear in atomic physics, acous- 
tics, quantum mechanics, solid state physics, aerodynamics, fluid mechanics, aeroa- 
coustics |1], |2], |8|, [9[, [15J, [16J. An important task of the theory of differen- 
tial equations is to find the minimal smoothness of coefficients, under which the 
equation admits a solution. Although there are some models that are closely re- 
lated to the differential operators with distributions in coefficients it is impossible 
to construct the theory of linear differential equations with distributional coeffi- 
cients, since the space of distributions T>'(R") is not an algebra with respect to 
the "pointwise" multiplication. This raises the basic question how to interpret the 
differential operators with distributions in coefficients. 

A lot of models are expected to be "selfadjoint" in the sense that appropriate 
operators, describing these models must be selfadjoint in some Hilbert spaces. Let 
a differential expression S correspond to such a model and let it contain distri- 
butions, supported by x — 0, in coefficients. In order to interpret the operator 
S we first construct a symmetric operator So by restricting S to the set of func- 
tions vanishing at the origin along with their derivatives. Then we consider all 
selfadjoint extensions of So and choose one of these extensions as a definition of 
S. This method goes back to the work of F. Berezin and L. Faddeev [6]. In some 
instances the set of all selfadjoint extensions of a symmetric operators is multi- 
parametric. Therefore the harder question comes: how to choose an extension that 
is best suited to our physical model. For some models the proper operator can 
not be chosen within the selfadjoint extensions theory, because the models con- 
tain hidden parameters. After replacing the singular coefficient with a sequence 
of short-range smooth coefficients, the operator obtained in the zero-range limit, 
as often happens, can depend on the type of regularization, i.e., the operator is 
governed by the shape of squeezed coefficients. This shape is a hidden parameter 
and plays a crucial role in the choice of a selfadjoint extension corresponding to 
the physical model under consideration. 

l 
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In [ 10], [11] the problem how to define the one-dimensional Schrodinger opera- 
tor with the ^'-potential, where S' is the first derivative of the Dirac delta-function, 
was considered . A natural approach to defining such a Hamiltonian is to ap- 
proximate 5' in 2?'(]R")-topology by regular potentials and then to investigate the 
corresponding family of regular Schrodinger operators. Therefore the authors con- 
sidered the family of Schrodinger operators on the line of the form 

approaching a formal Hamiltonian H a = —j^i + U(x) + kS'(x). Here e is a small 
positive parameter, Y 6 Cq°(— 1, 1), U is a real valued potential going to +oo as 
\x\ — > oo, and a is a real coupling constant. The map assigning a limit operator 
to each pair (a, Y) was constructed there. The choice of T-L(oi, Y) is de- 
termined by proximity of its energy levels and pure states to those for the Hamil- 
tonian with regularized potentials for small e. For almost all a the limit opera- 
tor is just the direct sum of the Schrodinger operators with the potential U on 
half-axes subject to the Dirichlet boundary condition at the origin (the nonresonant 
case). But for a belonging to the discrete resonant set which is the spectrum 
of the Sturm-Liouville problem — iv" + aWiv = on the interval ( — 1,1) subject 
to the boundary conditions w'( — l) = w'(l) = 0, the operator %(a,Y) acts via 
'H(a., x f)f = —f" + Uf on an appropriate set of functions satisfying the matching 
conditions /(+0) = 0yM/(-O) and %(a)/'(+0) = f'(-0) (the resonant case). 
Here 9^ (a) = w a {l) /w a { — 1), where w a is an eigenf unction corresponding to the 
eigenvalue a £ Hy. 

In | fl2) the results of [10], [11] were extended to the case of the fourth order 
ordinary differential operators. An attempt was made to define the formal dif- 
ferential operator A a = jj^ + U(x) + a.S"'(x). To approximate A a , the family 
A(«,Y) = Jjr + U(x) + £Y(f ) with domain 

WJ)) = {/€ W*(a,b): /(«)=/(«)=0, f(b)=f(b)=0} 

was considered. Here (a, b) is an interval of 1R containing the origin, U is a smooth 
real valued function on [a, b] and Y £ C^°( — 1,1). Asymptotic expansions for 
eigenvalues and eigenfunctions of ^4 £ (a,Y) were constructed, and therefore the 
limit operator A(oc, Y) was obtained. Upon constructing asymptotics two different 
cases are distinguished: the resonant case and the nonresonant one. In the resonant 
case a belongs to the discrete resonant set C 1R, which is the spectrum of the 
eigenvalue problem 

fll . iivW +ctYw = 0, ££(-1,1), 

1 ' ; ]y(-i) = w'"{-\) = o, w"{\) = w"'{i) = o. 

The limit operator was obtained under the additional assumptions 

(1.2) a. is a simple eigenvalue of the problem and w' a ( — l)w' a (l) ^ 0, 

where w K is an eigenf unction corresponding to the eigenvalue a. In this case 
*4(a,Y) acts via A.(&,Y)f = + Uf on an appropriate set of functions obey- 
ing the interface conditions f(0) = 0,/'(+0) = T (a)/'(-O) and &w(ct)f"(+0) = 
/"(— 0), where %(a) = w' a {l) / w' a ( — 1). In the nonresonant case, when a i 
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the limit operator is the direct sum of the Dirichlet operators on (a, 0) and (0, b) 
respectively. 

We extend the results of |[T2l to more a general perturbation of the fourth order 
differential operator, namely, we consider a formal differential expression 

d 4 

-j + U(x) + aS"'(x) + ps"(x) + yi5'(x) + l2 S(x). 

The investigation of the papers [10 1 — [12 1 is based on the asymptotic analysis. We 
will use techniques of [10[-[12[ to obtain an appropriate limit operator. 

1.1. Problem statement and main results. Let L stand for the differential expres- 
sion -j^g + U(x). As before If is a smooth real valued function on the interval 
[a, b] C K, containing the origin. Denote by Y £ the function 

Y £ (x) = ae^YlV 1 *) + jSe" 3 *^ 1 *) + 7 1 e" 2 Y 1 (£~ 1 x) + 7 2 £~ 1 Y 2 (£~ 1 x). 

Here Y, <3>, Y 1 ,Y 2 G Cg°(-1, 1), supp Y = [-1, 1], and a, 0, 71, 72 G 1R are arbitrary 
constants. Let us consider the eigenvalue problem 

(1.3) Ly £ + Y £ (x)y £ = A £ y £ , x G (a,b), y £ (a) = y' £ {a) = y £ {b) = y' £ (b) = 0. 

Note that the further analysis of the problem l|1.3t does not depend on the type 
of boundary conditions. Hence Dirichlet boundary conditions may be replaced by 
one of the possible combinations at the endpoints x = a and x = b of the following 
boundary conditions 

y(*o) = y'( x o) = 0, y{x ) = y"(x ) = 0, y"(xo) = y'"(x ) = 0. 

We associate with the problem dl.3t an operator 

d 4 

5 e (a, j S,7i / 72;Y,0 / Y 1 ,Y 2 ) = + U{x) + Y e (x), 

P(5 £ ( a , J 6,7 1 ,72;Y,* / Y 1 ,Y 2 )) = {/ e W 2 4 ( fl ,b): 

/(«)=/'(a) = 0, /(fc)=/ / (fc) = 0}. 

We denote it briefly by S £ . 

Note that for some Y,<E>,Yi,Y2 G Cg°(— 1,1) the function Y e converges in the 
sense of distributions as £ — >■ to the linear combination of the derivatives of the 
Dirac delta-function, which serves as a motivation for the choice of the singular 
perturbation Y £ . If therefore the operator S £ converges (in some sense) as £ — > to 
the limit operator, then it is natural to regard this limit as the interpretation of the 

fourth order differential operator ^ + U + ccS'" + f>&" + 71 J' + 72^. 

Our purpose is to investigate the asymptotic behavior of eigenvalues A £ and 
eigenfunctions y £ as £ — » 0. The perturbation Y £ consists of four terms each of 
which has different influence on A e and y e . It is of interest to know when each term 
starts to have effect in asymptotic expansions. Intuitively, we expect that the term 
approximating the third derivative of the Dirac delta-function has to be dominat- 
ing. We also wish to assign an operator to each collection (a, /3, 71, 72; Y, <I>, Yi, Y2) . 
We base the choice of the limit operator on the proximity of its eigenvalues and 
eigenfunctions to those of the operators S £ for sufficiently small £. 

The rest of the paper is organized as follows. Sec. 2 includes the description 
of the spectrum of the perturbed operators S £ . We show that all eigenvalues are 
continuous functions of £ and are bounded from above. Generally speaking, the 
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spectrum of this family is not bounded from below as e — > 0: for some Y and a 
there exists a finite number of eigenvalues converging to — oo as e — »■ 0. 

Then Sec. 3 presents the formal asymptotic expansions for the eigenvalues and 
eigenf unctions of S £ . The leading terms of asymptotic expansions and the limit 
operators are constructed in the section. We introduce a spectral characteristic of 
the shape Y, namely, the resonant set Y^, which is the spectrum of the eigenvalue 
problem In the case when a. does not belong to the resonant set, the limit 

operator is just the direct sum of the Dirichlet operators on (a, 0) and (0, V) respec- 
tively. In the resonant case, when a E Tpf, the limit operator <S a «(Y, <E>) acts via 
<S a ^(Y, <!>)/ = Lf on a set of functions obeying appropriate coupling conditions 
at the origin. 

The remainders of asymptotics for eigenvalues and eigenfunctions of S £ are 
constructed in Sec. 4, because we are in need of more precise asymptotics in order 
to prove the approximation theorems. In this section we also analyze the effect of 
each singular term. The justification and estimation of the range of validity for the 
approximations are presented in Sec. 5. 



in P'(1R) -topology. Set (f) k = (fc!)" 1 /"^ £*/(£) d.g. It is easy to prove that a func- 
tion / € C§°(-1, 1) is the <>»-like shape if and only if (/) ; - = for ; = 0, . . . , n - 1 
and (/)« = (—!)"( see [11 J for details). In what follows, we denote by M n the set 



For all e > the spectrum of S £ is real and discrete. Let {A|}^Lj be the eigenval- 
ues of S £ enumerated in increasing order taking multiplicity into account. Suppose 
that {y!.}j£L-i is the L^fl, b)-orthonormal system of eigenfunctions. 

Theorem 2.1. The eigenvalues Af. of the operator S £ are continuous functions of 
e £ (0,1). Moreover, all eigenvalues are bounded from above as £ — » 0. Let Y 
change sign and a \ be large enough; then the spectrum of S £ is unbounded from 
below as £ — >■ 0, in particular, < — ce -4 for some positive constant c. There is at 
most a finite number N~ of eigenvalues converging to — oo as £ — > 0. 

Proof. Let us consider the quadratic form 



that is equicontinuous on the set of functions u G T>(cj e ) fl {\\v\\ = 1} with respect 
to £. The minimax principle (5J p. 343] 

A| = inf sup q £ [v], 

fc * v£E k ,\\v\\=l 

yields continuity of eigenvalues with respect to £. Here runs over all A:-dimensional 
linear subspaces of H, and || • || denotes the 12(0, fc)-norm. 



2. Spectrum of S £ and auxiliary results 
An element / of C~(— 1, 1) is called the <!K") -like shape if 
g -(«+i)/( g -l x ) (jc) as £^0 





ueH = {fe W*(a,b): /(«)=/'(«) = (), /(&)=/'(&) = ()}, 
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Choose a subspace Et containing only elements vanishing in a neighborhood 
of the origin. Then we obtain 

A| < sup q E [v}. 

i>e£|, |M|=i 

For sufficiently small e the restriction of q E to £^ does not depend on e. This yields 
boundedness of the eigenvalues from above. 

Suppose Y changes sign. Let u £ C^{a, b) be a normalized function supported 
on an interval [ci,c 2 ], where Y takes negative values. Consider the sequence 
u £ (x) = £~ 1//2 m(£ _1 x), ||« £ || = 1, and assume that a > 0. From the minimax 
principle one can conclude that 

e 4 Af <e%[u E ] =£ 4 ^ C2E (|i4'| 2 + !J|m £ | 2 +Y £ | U£ | 2 ) dx = 

= £ 2 (| U "| 2 -«|Y|| U | 2 ) 

+ £ jf* (jSO + TieYx + 72 £ 2 Y 2 ) |m| 2 d£ + £ 4 jf* U(e£) |w| 2 d?. 
The first integral gives a negative number for cc > r, where 

r = £V(£)| 2 ^(/V(£)ll"| 2 ^) , 

while the other terms go to zero. Thus for £ sufficiently small the estimate \\ < 
— C£~ 4 holds with some positive c. The case a < may be handled in much the 
same way. 

Let Ng~ denote the number of negative eigenvalues of the operator S £ . Clearly, 
N~ < lim sup £ _j,Q N £ . It is well known |4J that the estimate for the number of 
negative eigenvalues 

N~<c + ci/' \x\ 3 \U(x)\ dx + c 2 [ \x\ 3 |Y £ (x)| dx, 

J a J a 

holds, where Cq, C\ and c 2 are positive constants. The function Y £ is supported on 

[ — £,£], thus 

j\x\ 3 \Y £ (x)\dx = J" |x| 3 (e _4 |aY(e _1 jt)| +£- 3 \p®(£- 1 x)\ + 

+ £- 2 |7iY 1 (£- 1 x)| +£- 1 | 72 Y 2 (£- 1 x)|) dx = \ tt \ J 1 |£| 3 |Y(£)|d£+ 

+ £ 1^ |?| 3 (|/5O(0| +£|?iY 1 (OI +e 2 |? 2 Y 2 (Ol) ^ < |«|c(T) + l 
for small £ > 0. From what has already been proved it follows that 

N~ < c + Cl (LT) + c 2 (Y)|a| 
for some positive constants Cq, Ci(U), c 2 (Y). □ 

Therefore the spectrum of S E consists of two parts: the set of eigenvalues tend- 
ing to —oo as £ — »■ 0, and the set of all bounded eigenvalues as £ — > 0. 
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3. ASYMPTOTICS OF EIGENVALUES AND EIGENFUNCTIONS OF S E AND THE LIMIT 

OPERATOR 

Fix an eigenvalue A| of the problem ( |1.3t with k > N~. We write it A e for short. 
Let y £ be the corresponding eigenfunction. The asymptotic expansions of A £ are 
represented by 



(3.1) 



A £ ~ A + eAj + £ A 2 + . . . , 



and we postulate two-scale expansions for the eigenfunction 

(3.2) y £ (x) ^v(x) + evi(x) + e 2 V2(x) + . . . for x G (a, — e) U (e, b), 

(3.3) y £ (x) ~ eiv(e~ 1 x) + £ 2 zui(e~ 1 x) + . . . for x € (— e,e). 

Here all functions v, are defined for x £ (fl, 0) U (0, i>), and if, are defined 
for £ 6 (—1, 1). Assume that w is different from zero. Series 113.21 1, d3.3b satisfy the 
coupling conditions 



(3.4) 



0, 



,0T 



0, j = 0,...,3, 



where by [/]x=« we denote the jump of / at a point a. 

We substitute series i3.1\ , §3.21 into the equation and the boundary conditions 
fll.3l l and derive 

(3.5) Lv = Av, v(a) = v'(a) = 0, v(b) = v'(b) = 0, 

(3.6) Lvi = Xv l +X x v, v t (a) = v[(a) = 0, p x (b) = v[(b) = 0, 

(3.7) Li; 2 = Au 2 + Ai^i + A 2 u, i» 2 («) = » 2 ( a ) = °> v 2 (b) = v' 2 (b) = 0, 

where all equations hold on («,0) U (0,b). We set £ = e" 1 *. After substituting 
d3.lt , ( 13.31 1 into the equation 11. 3\ , one obtains the following equations on ( — 1, 1) 



(3.8) 
(3.9) 
(3.10) 
(3.11) 



a,(4) + a Yz^ = 0, 

wf^ + aYw 1 = -p<&w, 

+ ctYii>2 = — fi&w\ — j^Yiw, 

.,(4) 



' + aYzt>3 = — /30ztf 2 — 71Y1W1 — 7 2 Y 2 zt>. 
Substituting Il3.2l l, 03 .3I > into the coupling conditions l|3.4|l , we can assert that 



v('\±e) + evP(±e) + e 2 v^(±e) + e 1 ^' »W(±1) + £^'^(±1) 



,2 „(/) 



2-/ „,(/), 
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for j = 0, ... ,3. We can now expand into the formal Taylor series about x = 
±0. Then we conclude that 



(3.12) 


v(-0) = 0, 


v (+0) = o, 


(3 13) 


w"(±l) = 0, 


w"'(±l) = 


ex 1 a\ 


V {—U) — W [ — 1 ), 


V I^tU ) — W [l ), 


(3.15) 


»l(-0j - v (-0) = w(-l), 


ZM+0) + v (+0) = w{l), 


(3.16) 


v[(-Q)-v"(-0) = w' l (-l), 


v[(+0)+v"(+0)=w' l (l), 


(3.17) 






v 2 (-0) 


-v' l {-0) + \v"(-0) = w l {-l), 


P2(+0)+oi(+0) + X(+°) = 


(3.18) 






v' 2 (-0) - 


-v l i{-Q) + \v'"{-Q) = w' 2 {-l), 


i4(+o) + i4'(+o) + !*/"(+o) = 


(3.19) 


w'({±\) = v"(±0), 


<(±1) = 0, 


(3.20) 


w' 2 '(±l) = v'{(±0)±v"'(±0) / 


<(±1) = p'"(±0), 


(3.21) 






a#(±l) = 




o#'(±l) = ?i"(±0)±i>( 4 )(±0). 



It follows that ?7 satisfies the equation and the boundary conditions (|3.5b , and fur- 
thermore v (0) = 0. The function w is a solution to the problem 

(3.22) a;( 4 '+aYa) = 0, fe(-l,l), a>"(±l) = 0, o/"(±l) = 0. 

Moreover these functions are related by the coupling conditions (|3.14b . The prob- 
lem d3-22|) is decisive in our next consideration, because it contains information 
about the singular perturbation. The first and primary question is whether there 
exists its nontrivial solution. 

3.1. Resonant set. The problem l|3.22b can be regarded as a spectral problem with 
the spectral parameter a. We note that in the generic case, the function Y is sign- 
changing. It is of interest to investigate spectral properties of this problem. We 
will also introduce the spectral characteristic of the shape Y. 

We introduce the operator 7% = fj^) JFi w i m the domain 

T>(Ty) = {fe Lf Y |(-l,l) | y-V (4) e Lf Y |(-l,l), /"(±i) = f"'(±l) = 0}. 

The problem l|3.22b is equivalent to the spectral equation T^w = —ctw. 

Theorem 3.1. Given Y 6 Q°(R), with suppY = [-1,1], the spectrum of the 
operator 7% is real and discrete. Suppose Y changes sign; then the spectrum of 7y 
has two accumulation points -co and +oo. 

Proof. Since the case where Y keeps sign is much simpler and can be handled 
within the standard Hilbert space theory, we assume that Y changes sign and 
apply the Krein space theory to investigate the spectrum of 7y- Let C be the 
weighted 1,2-space with the scalar product (f,g) = f \ \Y\fgdg, and let us de- 
fine by [f,g] = J_ 1 ~Yfgd£ the indefinite metric in C. Then the pair (£, [■, •]) is 
called a Krein space fl3j ch. 1]. 
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In this Krein space there exists the fundamental symmetry // = sgnY/ such 
that \f, g] = ( //, g) for all /, g G £. An operator T is J-selfadjoint if JT is selfadjoint 
in LL|(— 1, 1). An operator T is said to be J-nonnegative if [Tf,f] > for all / G 
V(T). 

For each Y G QJ°(— 1, 1) the operator /7% is selfadjoint, and so 7% is /-selfadjoint. 
Next, for all / G 2?(7y) one obtains 

[r Y /,/] = j\f®7<% = j\ I/" I 2 # > o- 

Thus 7% is /-nonnegative. Any /-selfadjoint and /-nonnegative operator with a 
nonempty resolvent set has real spectrum [3, p. 138]. Let us show that the resol- 
vent set of 7% is nonempty. The homogenous problem 

(3.23) g (4) +zTg = 0, £G(-1,1), g"(±l)=Q, g"'(±l) = 

has a trivial solution only. Indeed, each solution satisfies the equality 

£ i \g»\ 2 d£ + i£Y\g\ 2 dS = 0. 

Since Y is real-valued, it follows that g is a linear function. Obviously, only zero 
function can be a solution of d3.23l >. Hence the nonhomogeneous problem g^ + 
Vfg = h, g"(±l) = 0, g"'(±l) = admits a unique solution for arbitrary h G 
L 2 (-l,l) [13, p. 39]. Note that Y/ belongs to L 2 (-l,l) for each / G £, since 
IIT/1Il 2 (-1 1) — max R I 1 ?) 1 '' 2 ■ ||/|| £ . Then the equation Tyg + ig — f is equivalent 
to the nonhomogeneous problem g( 4 ) + zYg = Y/, g"(±l) = 0, g ;// (±l) = and 
admits a unique solution for each f G jC. Therefore —i belongs to the resolvent set. 
Since the resolvent set of 7% is nonempty, the spectrum of 7% is real. 

We shall prove that the resolvent R^(7y) of the operator 7y is compact. The 
operator K^(7y) acts from the space £ into V(Tk), and for each / G £ solves the 
equation 

g W - jffg = Y/, g G V{Tw). 

As far as the right-hand side Y/ belongs to L 2 (— 1,1), it follows that the solu- 
tion g is an element of W| ( — 1, 1). The space D(7y) is a Banach space with the 
graph norm. The sequence of continuous embeddings T>{ 7%) C Wo (—1,1) C 
— 1,1) C £ yields the compactness of the resolvent, since — 1, 1) C L 2 ( — 1, 1) 
is the compact embedding. As a consequence we have cr(Ty) = <7p(7r). 

Since Y changes sign, the spectrum <j(Ty) is unbounded in both directions [7J. 

□ 

We introduce the set T^g = {a G 1R : — a G c(7y) }/ which is the spectrum of the 
problem | |3.22|| . We call the resonant set of the shape Y. When a G Ly, suppose 
that \1.2) holds (the case of nondegenerate resonance). In this paper we assume that 
only the nondegenerate resonance is possible, namely, if a belongs to the resonant 
set, then both conditions l ll.2|l hold. 

3.2. The limit operator. Let us continue to construct the asymptotics. We distin- 
guish two different cases and start with the assumption a £ • Then the problem 
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113.221 admits a trivial solution w — only. That v'(0) = follows from the cou- 
pling conditions l(3.14|> . We conclude from B.51 that v is a solution to the problem 



(3.24) 



\Lv = \v, x G (a,0)U(0,Z>), 

ja(a) = p'(a) = 0, v(0) = v'(0) = 0, = p'(fc) = 0. 

Let us introduce the operators 

S_/ = L/, P(S_) = {/e W 2 4 (a,0): /(«)=/(«)=0 / /(0)=/'(0) = 0}, 

S+/=L/, D(S+) = {/eW 2 4 ((U): /(0)=/(0)=0, /(b) = /(&) = 0}. 

The operator S_ © S+ is associated with the problem ( |3.24t . Therefore in the non- 
resonant case, when a. £ Ey 7 we can define the limit operator as S_ © S+. 

Let us now suppose that ol belongs to the resonant set Ey- Recalling dl -2b , we 
deduce that the quotient 

<M a ) — — 

is well defined and does not depend on the choice of an eigenfunction. Clearly, 
w — cw K (£), where c is a constant. We conclude from J3.14I I that v'{— 0) = cu^(— 1), 
d'(+0) = cw^(l), hence that 

(3.25) z/(+0) - 0y(«)i/(-O) = 0, 

and also that c = ^p^ph ■ According to d3.9b , $3.19) the next term W\ of series i3.3i 
can be found by solving the problem 

w\ 4) + a x ¥w 1 = - 4^-$^, I G (-1, 1), 
(3 26) 1 fl ' ^ \ > i> 

w'{{-\) = o"(-0), <'(-*) = 0, w'{(l) = v"(+0), w'{'(l) = 0. 
Because a is an eigenvalue of J3.22I I, the problem admits a solution if and only if 

(3.27) M*K(+0) " ^'(-0) = ^'(-0) £ (^S)) 2 

To derive this solvability condition, we multiply the equation by w a and integrate 
by parts. Let us define a functional on Cg°(-1, 1) by &#[f] = j\<$>f 2 d£. Col- 
lecting i3.5) , d3.25|) and i3.27) we deduce that v must be an eigenfunction of the 
problem 



(3.28) 



'Lv = Au, x G (fl,0) U (0,b), 
z>(a) = o'(fl) = 0, v(b) = v'(b) = 0, 
p(0) = 0, i/(+0) - 0y(a)i/(-O) = 0, 
,0y(«)o"(+O) - o"(-0) - ^c[i^/<(-l)]^(-0) = 0. 

,4 

Consequently, the operator S(a, f>; Y, <J>) = ^4 + II (x) with the domain 

(3.29) 27(S(* / ^;T,*)) = {/e W 4 ((a,0) U (0,fc)) : /(a) =/'(a) = 0, 
/(6)=/(fe)=0, /(0)=0, /'(+O)-0 T (*)/'(-O) = O / 

(+0) - /" (-0) - ^ [ Wa /<(-l)]/'(-0) = 0} 
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is associated with the problem ||3.28|| . Combining resonant case and nonresonant 
one, gives us the limit operator 



S- © S+, mi E Y , 
S(a / j8;T,<E>) / a e E Y . 



Recall that we consider only those a from the resonant set, which satisfy assump- 
tions Oil. 



4. Asymptotic expansions of eigenvalues and eigenfunctions of S £ : 

CORRECTORS 

In order to justify the closeness of eigenvalues and eigenfunctions of operators 
<S £ and <S a/j g(Y, <&) we must derive next terms of series d3 . H > — 03.3b . Clearly, the 
construction of correctors depends on the multiplicity of A. Let A be a simple 
eigenvalue of 5 a ^(Y, O) with the eigenfunction v being normalized in L^ia, b). 

4.1. Asymptotics in the nonresonant case. In this subsection we assume that a 
does not belong to the resonant set Ey- Then a; = and c(5 a/ ^(Y, <!>)) = c(S_) U 
cr(S+). If A is a simple eigenvalue of <S a « (Y,$?), then A is a simple eigenvalue 
of S_ or S+. Without loss of generality we may assume A 6 c(S+), and thus i» 
vanishes on (a, 0). Employing | |3.9t , J3.19l > gives us the problem 

w[ i] + ct y ¥zv 1 = 0, £e(-l,l), 

<(-l) = <(-l) = 0, w'i{\) = v"(+0), <(1) = 0, 

which admits a unique solution, since a does not belong to the spectrum of ||3.22| |. 
In light of l!3.6l l, (|3.15l l the function v\ can be found by solving problems 

Lv\ = Xr>\, x G (a,0), 
fi(fl) = v[(a) — 0, 
[ Ul (-0) = 0, i/ 1 (-0) = a/ 1 (-l) / 

Lv\ = Ai»i +Aii7, x 6 (0,fc), 
Dl (+0) = -u'(+0), v[(+0) = - v"(+0), 
[ Vl (b) =v[(b) =0. 

on (a, 0) and (0, b) respectively. Of course, the first of these problems has a unique 
solution, since A i cr(S-). Note that in the generic case the second problem has 
no solution. But we can ensure the existence of a solution by choosing the free 
parameter X\. Indeed, applying the Fredholm alternative we conclude that the 
second problem ||4.1| | admits a solution if and only if 

Ax = v"(+0)(v"(+0) - w[(l)) - z/t+OK't+O). 

To derive this we multiply the equation by the eigenfunction and integrate by 
parts. The last equality is simultaneously a formula for the corrector Ai in the 
asymptotic expansions of the eigenvalue. Clearly, the solution V\ is defined up to 

the term cv. To fix it we subordinate the solution to the condition J Q vv\ dx = 0. 



(4.1) 
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(4.2) 



Combining j3.10\ with J3.20[) and recalling a £ Ey, we deduce the problem 

w { 2 ] + a x ¥w 2 = -p><$>w lr I e (-1,1), 
w' 2 \-\) = ^'(-O), <(-l) = 0, 
< (1) = »i'(+0) + i/"(+0), <(1) = r/"(+0), 
which gives us the corrector w 2 . We employ < |3.7> , d3.17|) . <3.18fr to find 

' Lv 2 = Xv 2 + AiT7i, ie (d,0), 
J7 2 (fl)=0, i> 2 (a) = 0, 
p 2 (-0) = ipi(1) + o' 1 (-0), 
Ui(-0) = B)' 2 (-l) + r»i'(-0), 

' Lj?2 = hv 2 + A-jOj + XqV , x 6 (0, b), 
i> 2 (+0) = ii^(l) - v[{+0) - \v"{+0), 
v' 2 (+0) = w' 2 (l) - v'{{+0) - \v"'{+0), 
{v 2 (b) = v' 2 (b) = 0. 

As before we deduce that the first of these problems has a unique solution, and 
the second one admits a solution if and only if 

A 2 = t/"(+0) (wi(l) - Di(+0) - \v"{+0)) - v"(+0) (w' 2 (l) - v'{(+0) - \v"'{+0)). 
For the sake of definiteness, the solution is subject to the additional condition 
Jq vv 2 dx = 0. By using ( |3.1H . I|3.21t one obtains 

xv^ + ocYw 3 = -fi<S>iv 2 - 7j Y x wi, £ G (—1,1), 

«#(-!) = ^'(-0) - <(-0), <(-l) = <(-0), 

^'(1) =4'(+0)+z;i"(+0) + iz;( 4 )(+0) / <(1) =<(+0) + i;( 4 )(+0). 

Reasoning as before, from this problem we get zt>3. 
Let us introduce the notations 



(4.3) 



Ve(«) 



A e = A + eAi + e 2 A 2; 

i?(x) + eci(x) + e 2 i> 2 (x), x £ (a, — e) U (e, fr), 

£ 2 wi(e _1 x) + £ 3 zf2( e_ lx ) + £ i w 3 {t~ 1 x), x g (— e,e) 



for the constructed approximations of eigenvalues and eigenfunctions. 

4.2. Asymptotics under resonance. Now we assume that a belongs to the reso- 
nant set Ey and that A is an eigenvalue of the operator S(a, /5, Y, <J>). Let w a be an 
eigenf unction of l|3.22b such that zi^( — 1) = 1. Clearly 0y(«) = iy a(l)- 

Since | |3.27| | holds, the problem d3.26D admits a solution. This solution can be 
represented as W\ = w\ + C\Wm, the function w\ being a partial solution of the 

problem fixed by the condition 1) = 0, and the constant c\ is to be chosen 

later. 



We next construct the corrector V\. The function v\ satisfies the equation 1 13.6 
outside the origin and (T3l9t yields 



(4.4) 



v[(+0) - e Y (cc)v[(-0) = Gi, 
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where G r = w[(l) - 9w{oc)w[(-l) - v"(+0) - &y{oc)v"(-0). Although w x is not 
uniquely chosen, the constant G\ is well defined. In fact, 

w[(l) - M«K(-i) = (^f(i) " + 

+ Cl «(i)-M«K(-i)) = ^(i)- 

From (|3. 10b and d3-20b it follows that the corrector Wi must solve the problem 

+ of¥w 2 = -$®Wi - 7 1 Y 1 w, £ G (-1, 1), 

(4.5) ^'(-1) = ^(-o) - v'"{-0), <(-l) = o"'(-0), 
h#(1) = v'{(+0) + z/"(+0), <(1) = i/"(+0). 

From the first condition in <|3.16b we deduce c\ = v[(— 0) — v"(—0). Set i?Y-[f] — 
/* x Yjfd£ for / e Cg°(-l,l) and z = 1,2. Thus the solvability condition of the 
above problem can be written as 

(4.6) M«K(+0) - n?(-0) - 0*rKK(-O) = Hi 
which is due to the Fredholm alternative. Here 

Hi = - m*)K'(+o) - + i>"'(-o)+ 

+ ^[^w*^\ -^KK't-oD+Ti^Kjs't-o). 

From J3- 15|> we have »i(±0) = 0)w a (±l) =F »'(±0). Combining these identi- 
ties along with l|3.6|l , l|4.4|l and l|4.6t we conclude that V\ solves the problem 
(4.7) 

{Lvx = Xvi + Ait?, x e (fl,0) U (0,b), 
v x (a) = v[{a) = 0, 17! (b) = ^(fc) = 0, 
*i(-0) = p'(-0K(-1) + v'(-0), v 1 (+0) = v'(-0)rv K (l) - v'(+0), 
v[(+0) - 9 w {oi)v[(-0) = Gi, 9 w {a)v'{(+0) - v'{{-0) - ^[w x ]v[(-0) = H v 

The free parameter Aj in the right-hand side of equation J4.1 lb enables us to solve 
the problem. In view of Fredholm's alternative, J4-1 1|) admits a solution if and only 
if 

Ai = Hii7'(-0) - G^"(+0) - (i/(-0K(-l) + v'(-0))v"'(-0)+ 

+ (v'(-0)w a (l) -v'(+0))v"'(+0). 
For the sake of definiteness, the solution is subject to the additional condition 

la VV l dx = °- 

Given v\, we may compute the constant c\. A trivial verification shows that the 
second condition in B.16|l holds. 

The condition < 14.6l l enables one to solve the problem | |4.5l l. A solution of this 

problem has the form w 2 = + c 2 Wol , where w\ solves 1)4.51 and satisfies ( — 1 ) = 
0. The constant c 2 will be chosen later. 

Following as before we shall similarly find correctors ZV3, v 2 and A2. The func- 
tion v 2 satisfies the equation and boundary conditions i3.7) , and 

(4.8) v' 2 (+0) - 9 w (a)v' 2 (-0) = G 2 
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by&M- HereG 2 = ^(1) - T («) «(-0) - ^'"(-0)) - ^'(+0) - ^'"(+0). 

Next we employ l[3.11t and d3.21|) to obtain the problem 

(4.9) 

tflg + afwj, = —p<&w 2 — 71Y1W1 - 72Y2W, £ e (-1, 1), 

w'i(-\) = v'l(-O) - v'{'(-0) + ^( 4 )(-0), w'i'{-l) = if(-0) - v^(-O), 

a#(l) = v'i(+Q) +v'{' (+0) + ^ (4) (+0), a#'(l) = *f (+0) W 4) (+0). 

On applying 0.18|l , one obtains c 2 = v' 2 { — 0) — i?"(— 0) + lv"'(— 0). Therefore we 
may write the solvability condition of this problem in the form 

(4.10) m«K(+o) - Pz(-0) - ^[w«]^(-0) = H 2 

with 

H 2 = v'{'(+0){w K (l) - 0*(*)} - v'{'(-0){w a (-l) + 1} + 

+ v W(+0){w x (l) - \6^{a)} + vW{-0){w K (-l) + + 
+ p(&r[ y JwZw eL ] - K(-0) - l/'(-0)}i%K]) + 7i# Yl [wiw a ]+72%[ww a ]. 
From | |3.17| | we find »2(=t0) = E^' wnere 

Ef = W J(±1) + (^(-0) - i;"(-0)) H ; a (±l) =f ^(±0) - ^"(±0). 
In view of $37), (SJ) and i HJO) it follows that c>2 is a solution to the problem 

' Lv 2 = A»2 + + A 2 w, x e (a, 0) U (0, fr), 
i> 2 (a) = w 2 («) = 0, v 2 (b) = v' 2 (b) = 0, 

= F 2 ' = F 2 > 4(+0) - M«K 2 (-°) = G 2/ 

.M*K(+0) - i#(-0) - ^[wJo^-O) = H 2 . 
The problem admits a solution if and only if 

A 2 = H 2 v'(-0) - G 2 v"(+0) - F 2 -v"'(-0)+F+v"'(+0). 
This solution is defined up to the term cv. To eliminate this ambiguity we addi- 
tionally demand that the condition J vv 2 dx = holds. 

Summing up, one obtains the following approximations for the eigenvalue and 
eigenfunction of the perturbed problem in the resonant case: 
(4.12) 

A e = A + eAj + e 2 A 2/ 

v(x) + £V\(x) + e 2 v 2 (x), x 6 (a, — e) U (e, b), 

erw(£ _1 x) + £ 2 zvi(e~ 1 x) + e 3 a; 2 (£ _1 x) + e^w^e^x), xe (—e,e). 

Here wj, is an arbitrary solution of l|4.9t . The choice of C3 in the representation 
W3 = W3 + c^Wu is not important since we do not look for the corrector C3. 

5. Justification of asymptotic expansions 

As shown in Theorem l2.ll for every regularization Y E (x) there is at most a finite 
number N~ > of eigenvalues A£, converging to — 00 as e — > 0. Other eigenvalues 
remain bounded as e — » 0. We shall show that these eigenvalues converge to the 
eigenvalues of <S a «(Y, O). 



(4.11) 



Y £ (x) = 
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5.1. Convergence theorem. Let {A £ } £e j be a sequence of eigenvalues of S £ and 
assume that {y e } £G x is a sequence of the corresponding Li(a, fc)-normalized eigen- 
functions. Here X is an infinite subset of (0,1) for which is an accumulation 
point. 

Theorem 5.1. If A £ — > A and y £ — > v in Li{a, b) weakly as I 3 e — >■ 0, then A is an 
eigenvalue ofS X/ p(Y, O) with the corresponding eigenfunction v. Furthermore, y £ 
converges tov in L2 (fl, fo) . 

We have divided the proof into a sequence of lemmas. To start with, let us 
describe the behavior of y £ outside the £-neighborhood of the origin. 

Lemma 5.2. Under the assumptions of Theorem 15.11 for every positive 7 the se- 
quence y £ tends to v weakly in the topology of W| ((a, b) \ ( — 7, 7) ) and strongly 
in the topology ofC 3 ([a, b] \ (—7,7))- Furthermore, v solves the equation 

(5.1) Lv = Av, x G («,0)U(0 / fe). 

Proof. Throughout the proof Gj denotes the set of test functions cp 6 Qj°(a, fc) such 
that <p{x) = for x £ (—7, 7). From the equation d 1 -3b for all cp G G-y and £ < 7 we 
deduce 

rb rb 

(5.2) / Ly £ cpdx = A £ / y c cpdx, 

J a J a 

since supp Y E C (—7,7). The right hand side of l|5.2b has a limit as 1 3 e — > 
by assumption, thus the integral on the left hand side converges for all cp 6 G-y ■ It 
follows that y E — >■ i> in W| ( {a, b) \ ( — 7, 7) ) weakly and thus 

Lvcpdx = \ vcpdx, cp £ Gj- 



From this identity it may be concluded that v solves d5.lt on (a, b) \ (—7, 7), and 
so on (a, 0) and (0, b), since 7 is an arbitrary constant. Applying the imbed- 
ding theorem yields convergence of y £ in C 3 ((a, b) \ (—7,7)), which completes 
the proof. □ 

We proceed to investigate the behavior of y £ along with its derivatives at the 
points x = — e and x — e. 

Lemma 5.3. Let A £ — > A and y £ — > v in Li{a, b) weakly as X 3 e -3- 0. Then for 
k = 0, . . . , 3 the sequences yi' (±e) converge to (±0) as I 9 £ — > 0. 

Proof. Let £fc be C™((«, b) \ 0)-functions such that £/ c (x) = for x < and ^(x) = 

k u 

jy for x 6 (0, !)• Denote by X( E ,oo) t ne characteristic function of (e, 00) and set 

£*(*) = X(e/») (*)£*(*)• We note that J^£f(x) = for x e (e,|). Multiplying 
both equalities fll.3|l , j5.1|l by £q and integrating by parts yield 

yi"(e) = - y"'Co ^ + [ vMv - a e ) dx, 

b rb 



J b v"'C dx + £ vCo(V-\)dx. 
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The right hand sides of the equalities have the same limit as X 3 e — > in view of 
Lemma I5~2l and so y"'(e) — > i/"(+0). Applying the function £q( — x) similar to the 
above implies y"'(— e) — > v"'(—0). We have proved the Lemma for k — 3. 

The case k = 2 can be handled in much the same way, the only difference being 
in applying the function £| . Multiplying jl.3i , d5.H l by £j and integrating by parts, 
we derive 

,-b r b 



yi'(e) = y?'(e> " J b _ v"£l dx + jf y £ £i(A £ - V) dx, 
v"(e) = v"'(e)e - jl v"C{ dx + J* v& (A - V) dx. 



2 

Again employing Lemma [5721 we deduce u" (e) — > v"(+0) as X 3 e — >■ 0. Using 
£f(— x) instead of £? (x) yields y" (— e) — > v"{— 0). The rest of the proof runs as 
before. □ 

We denote by g\, gi solutions of the following problems 

4 4) + aT(Og/c = 0, £e(-l,l), 

gk(-l) = S u , g' k (-l) = 8 Z/h Si'(-1) = 0, *f(-l)=0, 
where <5, y is the Kronecker symbol. Let g £ solve the Cauchy problem on [—1,1] 

(5 4) 2 W + (aY(?) + e2 7iYi(^))g = (e - V(-e)si(£) + y£(-e)#(£)), 

g(-l)=0, g / (-l) = 0, *"(-!) = *"H>), g"'(-l) = 0, 
and let g| be the solutions of the problems 

gf ] + («Y(£) +£ 2 7iYi(^))ft = 0, £ e (-1,1), 

^(-1)=^, ^'(-1)=0, 0. 

The task is now to describe the behavior of the eigenfunction y £ in the £-neighborhood 
of the origin. 

Lemma 5.4. If A £ — >■ A andy £ — > y in L,2(a,b) weakly as I 3 e — > 0, then 
(5.5) |ye(-e)| + |y e (e)| <ce, 

and moreover, 



(5.6) 



->0. 

C3([-l,l]) 



Proof. First let us prove that the sequence w £ (£) = £ 1 y £ (e^) — £ 1 y £ (— e)gi(£) 
y E (-e)g2(f) tends to in L 2 (-l,l) as X 3 £ -> 0. Set W £ (£) = £ _1 y £ (££) 
£ _1 y £ (— e)h\(£,) — y £ (— e)/z|(£) with /z| being solutions to problems 

4 4) + (aY(£) + ep*(e))fcjfc = 0, £ 6 (-1, 1), 
fc*(-l)=*i* ^(-1)=%, ftt(-l) = 0, <(-l) = 0. 
By construction £Y £ solves the problem 

f U ( 4 ) + («Y(0 + e/5cl>(0)M = /z £ (a £e (-1,1), 
ltt(-l)=0, U '(-l) = 0, «"(-!) = £y £ '(-£), u"'(-l)=£ 2 y £ "( 
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where h e (g) = — ey £ (e£)(7iYi(£) + £72Y 2 (£) + £ 2 U(e£) — £ 2 A E ). For every func- 
tion h e 6 1,1) the solution W £ of the above problem is unique, belongs to 
W 2 ( — 1, 1), and satisfies the estimate 

ll^llw 2 4 (-i,i) ^ c INIl 2 (-u) 

with constant c being independent of £. Employing the inequality 

e f_ x y 2 M)d^ = f_y\{x)dx < \\ye\\l z{atb) = h 

one concludes that £ 1/2 \\y E (e^) ||l 2 (-i,i) < c , hence that ||?i £ ||l 2 (-i,i) — C£l/2 , an d 
finally that ||2^e|| vv 4 (— 1 1) ~~ ^ as ^ 3 e ~~ ^ 0- From this it follows that 

(5.7) £ -V(e) - e-VC-OMCl) " ^("^1(1) -> 0, 

(5.8) ^(e) -e-V(-e)(fe!)'(l) -^(-^(fcD'W ^ 0' 

(5.9) ^'(e) - £ -Ve(-e)(M)"(l) -^(-£)(^I)"(1) 0, 

(5.10) ^"(e) -e-Ve(-c)(M)'"(i) -ye(-£)(^)'"(i) -> o, 

By <5lfl- (51C)) we deduce 

y 8 (-e)(Af) (fc) ( 1 ) = °( e ) as ^9^0 for Jfc = 1,2,3. 

If at least one of the values {h\)"{\) or {h\)'"{\) is nonzero, then y £ ( — e) = O(e) 
as X 9 e -> 0. Suppose, to the contrary, that = {h\ )'"{\) = 0. From (53) 

it follows that is an eigenfunction of the problem J3.22D corresponding to the 
eigenvalue a £ Ey- By construction (— 1) = 0, contrary to (1.2) . The proof of 
(53) is complete by using (5J7). 

Finally since — ^|||l 2 (-i,i) converges to as X 9 e — > for A; = 1,2, then u £ 
converges to in L2(— 1, 1) as X 9 e — >■ 0. In fact, 

Il**ellta(-1A) - ll W ellL 2 (-l,l)+ £ ~V(- £ )ll /l l-^llk 2 (-l,l)+ 
as e — >■ 0. 

We introduce the function W £ (f) = £~ 2 y £ (££) — £~ 2 y £ ( — e)gf(£) — £ _1 y £ ( — £)g 2 (£) — 
g £ (f )■ Analysis similar to the above implies that W £ solves the problem 

M ( 4 ) + («Y(0+£ 2 7lYi(?))"=/e(a f e (-1,1), 
u(-l) = 0, u'(-l)=0, it" = y'/(-e) - i/'(-0), u"'(-l) = eyf (-e), 

with/ e (£) = — /3<J>w E (£) — ey £ (e^)(7 2 Y 2 + eli(e£) — eA E ) and satisfies the estimate 

||W £ || W 4 ( _ U) < c (||/ 8 || £2( _ u) + W:{-e) - v"(-0)\) 

with constant c being independent of e. According to Lemma 15.31 the right hand 
side of this estimate converges to zero. Thus W £ tends to in W| ( — 1, 1 ) as X 9 £ — > 
0. The embedding W 2 4 (-l, 1) C C 3 ([-l, 1]) establishes the convergence (52). □ 
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Letting l9£->0 we conclude that 

(5.11) e- 2 (y £ (e) - y £ (-e)g\(l) - ey £ (-e)£(l)) - g e (l) -> 0, 

(5.12) e-\y' £ (e) - ^(-e) (^)'(l) -y' e (-e)(g e 2 )'(l)) - g' e (l) 0, 

(5.13) y>> { e)- £ -iy £ (-e)(giy>(l)- £ ^yi(- £ )( g2 y>(l)- g >>(l) _> 0, 

(5.14) £ y £ "( £ ) - e -2y £ (_ e ) (gl )'»(i) -e-V e (-e)C&)'"(l) -^"(1) 0, 
in light of Lemma IBT41 for £ = 1. 

Lemma 5.5. Assume mat A £ — > A and y e — > y in £2(0, b) weakly as I 9 e — > 0. 
Then y £ — >■ y in £2(0, b) as 1 3 e — >• 0. 

Proof. First we show that y £ is bounded on [a, b] uniformly with respect to e. Ap- 
plying Il5.6b and Lemma 15.31 we see at once that the sequence y s is uniformly 
bounded on [—£,£] 

(5.15) |y e (^)l <ce 2 + |y e (- e M(?)|+e^(-£)gI(?)|+ £ 2 |^(?)| <c :£ . 

Set Q £ = (a, b) \ (—e,e). Multiplying the equation ( |1.3> by the function ^n E ye and 
integrating by parts give 

/ V?dx = J {X e -U) V ldx-y' e "{-z)y e {-e)+ 

y £ "( £ )y £ ( £ )+^(- £ )y £ (- £ )-y £ '( £ )^( £ ). 

All terms on the right-hand side are uniformly bounded with respect to e. Thus 
the sequence y e is bounded in W? (fi £ ), and so in C 1 (Q £ ). On account of the above 
conclusion combining with 05.151 1, we deduce that max |y £ (x)| < c with constant 

c being independent of e. 

Fix 7 > 0. According to Lemma [5.21 the difference y £ — y has the L2(Q7)-norm 
less than 7, provided e is small enough. Then 

he ~ y\\i 2 ( a ,b) < he - ylk 2 (n 7 ) + Wye - ylli«(- w ) ^ 

< 7(1 + 2 max \ye(x)-y{x) |) < C7 
xe(«,fe) 

with constant C being independent of £. Recall that 7 may be made arbitrary small, 
and the proof is complete. □ 

Proof of Theorem |5Jl We conclude from Lemmas 15.2115.51 that v is a solution of the 
equation 

Lv = Xv, x e («,0) U (0,&), 

satisfies the boundary conditions tf(a) = v'(a) = 0, = v'(b) = 0, and 

II 17 IIl 2 (<! 6) = Furthermore, v(0) = according to Il5.5b , and we are left with 
the task of showing that v satisfies appropriate coupling conditions at the origin. 

Again applying | |5.5l l, we deduce that £ _1 y £ (— e) has a limit as X 3 £ — > 0, which 
will be denoted by s, and y' £ (—e) — > v'(— 0) as shown in Lemma 15.31 Therefore 
<? e = £ -1 ye(-£)Si + y' £ {-£)g\ converges to q = sg x + v'{-0)g 2 in C 3 ([-l,l]) as 
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X 3 £ — > 0. From l|5.13l l- (|5.14l l it may be concluded that sequences e 1 £7" (1 ) and 
e _1 ^"(l) are bounded as X 3 e — >■ 0, hence that 

(5.16) <?"(!) = 0, <j"'(l)=0. 
Combining d5.3l l with H5.12| | yields 

(5.17) q'(-l) = v'(-0), q'(l) = v'(+0). 
We see at once that q is a solution of the Cauchy problem 



(5.18) 



?(-!)= s, <?'(-!) = t/(-0), «7"(-l)=0, <H-1)=0. 



Coupling conditions of the limit problem depend on whether the problem l|5.18t 
admits a nontrivial solution. Let us suppose for the moment that the problem 
115.181 has a trivial solution £7 = only. Then ([5.171 implies the coupling condition 

o'(0) =0. 

Next assume that J5. 1 8b has a nontrivial solution. By 1)5.161 q is an eigenfunction 
of the problem d3.22|) and a belongs to the resonant set Tif. In view of l|5.17t we 
have 

q l (l)v'{-0)-q'{-l)v'{+0)=0, 

which is equivalent to v'(+0) — 6y{ix.)v'(—0) = 0. 
For every e > the Lagrange identity holds 

Dividing the above identity by £ and letting X 3 £ — » 0, we derive 

(5.19) z'"(l)<7(l) + (p"(+0) - z"(l))^(l) = 0. 

in light of J5.13|) , J5 . 14b , where z solves the problem 

fzw + «t(?)z = e e (-1, 1), 

\z(-l)=0, z'(-l)=0, z"(-l) = v"(-0), z"'(-l)=0. 
Taking into account (|5.191 , the Lagrange identity for z may be written as 

z'"(lMl)-z"(lV(l) + ^(-0) (? '(-l) = (z"'(lMl) + ^"(+0)-z"(l))^(l))- 
-v"(+0)q'(l)+v"(-0)q'(-l) = -v"(+0)q'(l) + v"(-0)q'(-l) = -^[q]. 

Dividing the last equality by q' ( — 1 ) and recalling l|5.17t gives 

M«K (+0) - v"(-0) - [w u /w' a (-l)]v'(-0) = 0. 

Thus v is an eigenfunction of S a ^ ( Y, <I>) corresponding to the eigenvalue A. □ 

Theorem 15. II allows one to justify the choice of 5 a;j g(Y, <J>). 

Theorem 5.6. Suppose that the eigenvalue A £ ofS £ is bounded from below. Then 
A £ has a finite limit as X 3 £ -3- and this limit is a point of the spectrum of 
<Sfl, ^ For each simple eigenvalue A of S K «0?, <3>) there exist exactly one 

eigenvalue A £ ofS E converging to A as £ -3 0. 
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Proof. Suppose to start with that 

H* = !im E _> A' c < lirrig^o^l = Y* ■ 
The constants \i* are finite since A£ is a bounded function. Recall that Ajr is 
a continuous function of £ G (0,1)- Then for each A G there exists a 

subsequence of eigenvalues A £ , e G I, converging to A. The sequence {ye} £ ei of 
the corresponding normalized eigenfunctions contains a weakly convergent sub- 
sequence. By Theorem l5.11 A is an eigenvalue of S a/j g(Y, <!>). Therefore the interval 
belongs to the spectrum c(S ci r (Y, <J>) ), a contradiction. 
We now turn to the second part of the theorem. Let us assume that Ajr — >■ A 
and Ajr +1 — > A for some k. Then there exist two sequences {yl} E el and {y^+ileeZ 
of eigenfunctions, which converge in L 2 («,£>) to vectors of the form e l ^v. This 
contradicts the fact that y| and y| +1 are orthogonal in h%{a, b) for all £ G X. □ 

5.2. Approximation theorem. We proceed to show that each point of cr(S a ^ (Y, <5) ) 
is a limit of the eigenvalues of S £ . 

Let B be a self -adjoint operator in a Hilbert space H with domain T>(B). A pair 
(ff, m) G IR x T>(B) with ||w||h = 1 is called a quasimode of the operator B with an 
accuracy p > if ||Bm — f/«||H < i°- 

Lemma 5.7. Suppose that the spectrum of B is discrete and simple. If (pi, u) is a 
quasimode of B with accuracy p > 0, then the interval [pi — p, pi + p] contains an 
eigenvalue A of B. Furthermore, if the segment [pi — t, pi + t] contains only this 
eigenvalue of B, then \\u — v\\h < 1x~^p, where v is a normalized eigenfunction 
of B for the eigenvalue A. fl4) 

Let us construct the quasimodes of S E . Suppose A is a simple eigenvalue of the 
operator S a ^(Y, <I>) with the eigenfunction v such that ||c|| = 1. Here and subse- 
quently, || ■ || stands for the norm in L 2 (fl, b). For each A and v we have obtained 
the formal asymptotic approximations A E , Y e defined by either <|4.3b or (|4.12|l de- 
pending on a and Y. In further computation we do not distinguish the resonant 
and non-resonant cases. By construction we have 

(5.20) 

LY £ - A £ Y £ = £ 3 #i (e, x), x G (a, -e) U (e, b), 

(L + oaT^(er l x) + Be^Q^x) + 7 1 £~ 2 Y 1 (£~ 1 x) + 7 2 £~ 1 Y 2 (£~ 1 x)) Y E - 

— A £ Y E = £R 2 (£,x), x G (—£,£). 

The function Y E does not belong to the domain of S £ , since it has jump disconti- 
nuities at the points ±£. Indeed, 

[Ye]x=±e = e 3 4(e), [Yi] x=±e = e 3 rf(e), 

[Y':\ X=±E = e\±(e), [Y' t "] x=±E = e 2 rf(e). 

Here all the functions R,, are uniformly bounded with respect to their argu- 
ments. We can construct a function £ £ with the following properties 

o £ £ is a smooth function outside the points x = ±£ and differs from zero 
only for £ < |x| < 1; 

o [C e ]*=± £ = -er±(e), [C £ } x =±e = -of (e), [C £ 'h=±e = -erf(e) and [£"]*=± e = 
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o max (!£«(*) I + \C e (x)\ + + \C"(x)\ + K""(*)l) < c with constant 

£<|*|<1 

c being independent of £, 
which eliminates discontinuity. In fact, Y £ + £ 2 £ £ is a function from C 3 (a, fc) and 

belon gs to 2?(«S e ). We set ^ = \\Y e + e 2 £e|| _1 ( Y e + £l Ce) and substitute J, into 
1 15.201 1 instead of Y £ . Then the orders of smallness of right-hand sides in H5.20I I 
do not change since ||Y £ + £ 2 £ £ | — > 1 as £ — > 0. Therefore the pair (A e ,y e ) is a 
quasimode of S £ with accuracy £. 

Theorem 5.8. Given (a,jS,7i,72;Y,4>,Yi,Y 2 ) G M 4 x (Cg°(-1, l)) 4 , suppose thaf 
A is a simple eigenvalue ofS a p (Y, <J>) with fhe normalized eigenfunction v. Then 
there exists a simple eigenvalue A| of <S £ with a corresponding normalized eigen- 
function j/J such that 



(5.21) 



A 



with constants C\, Ci being independent ofe. 

Proof. Let (A £ , y e ) be a quasimode of S £ corresponding to the limit eigenvalue A 
and the eigenfunction v. According to Lemma 15.71 there exists an eigenvalue A^ 

such that A e - — A £ | < C\£, from which the first inequality in (|5.21t follows. In view 
of Theorem [52] the index ;' is independent of £. If t is less than the distance from A 
to the rest of the spectrum of S K> p(Y, O), then the interval [A — t, A + t] contains 
the eigenvalue AJ only, provided £ is small enough. Applying again Lemma 15.71 

yields yj — y e < 1t~^c\z, from which the second inequality in J5.21|l immedi- 
ately follows. □ 



References 

[1] S. Albeverio, F. Gesztesy, R. Hoegh-Krohn, H. Holden, Solvable models in quantum mechanics. 

With an appendix by Pavel Exner, RI: AMS Chelsea Publishing, 2005. 
[2] S. Albeverio, P. Kurasov, Singular perturbations of differential operators and solvable Schrodinger 

type operators, Cambridge, Univ. Press, 2000. 
[3] T. Azizov, I. Iokhvidov, Linear operators in space with an indefinite metric, Pure and Applied 

Mathematics, Chichester, 1989. 
[4] M. Birman, M. Solomyak, Estimates for the number of negative eigenvalues of the Schrodinger operator 

and its generalization, Adv. Sov. Math. 7 (1991), 1-55. 
[5] F. Berezin, M. Shubin, The Schrodinger equation, Kluwer Academic Publishers, 1991. 
[6] F. Berezin, L. Faddeev, A remark on Schrddinger's equation with a singular potential, Sov. Math. Dokl. 

2 (1961), 372-375. 

[7] B. Curgus, H. Langer, A Krein space approach to symmetric ordinary differential operators with an indef- 
inite weigth function, J. Diff. Eq. 79 (1989), 31-61. 

[8] Yu. Demkov, V. Ostrovskii, Zero-range potentials and their applications in atomic physics, 
Leningrad, Leningrad Univ. Press, 1975. 

[9] N. Goloscshapova, L. Oridoroga, 4-th order differential operator with local point interactions, Ukr. 
Math. Bull. 4 (2007), 355-369. 
[10] Yu. Golovaty, S. Man'ko, Schrodinger operator with 8' -potential , Dopov. Nats. Akad. Nauk Ukr., Mat. 

Pryr. Tekh. Nauky. 5 (2009), 16-21. 
[11] Yu. Golovaty, S. Man'ko, Solvable models for the Schrodinger operators with S'-like potentials, Ukr. Mat. 
Visn. 6 (2009), 173-206. 

[12] S. Man'ko, Fourth order differential operators with distributions in coefficients, Visn. L'viv. Univ., Ser. 

Mekh.-Mat. 72 (2010) (accepted for publication). 
[13] M. Najmark, Linear differential operators, Moscow, Nauka, 1969. 



SPECTRAL PROPERTIES OF THE FOURTH ORDER DIFFERENTIAL OPERATOR 



21 



[14] M. Vishik, L. Lyusternik, Regular degeneration and boundary layer for linear differential equations with 

small parameter, Usp. Mat. Nauk 12 (1957), 3-122. 
[15] A. Yablonski, Differential equations with generalized coefficients, Nonlinear Anal., Theory Methods 

Appl. 63 (2005), 171-197. 

[16] S. Zavalishchin, A. Sesekin, Dynamic impulse systems: theory and applications, Dordrecht, 
Kluwer Academic Publishers Group, 1997. 



Department of Mechanics and Mathematics, Ivan Franko National University of 
Lviv, 1 Universytetskastr., 79000 Lviv, Ukraine 
E-mail address: s_maiikoafranko.lviv.ua 



